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ABSTRACT 

High redshift galaxies permit the study of the formation and evolution of X-ray binary populations on cos- 
mological timescales, probing a wide range of metallicities and star-formation rates. In this paper, we present 
results from a large scale population synthesis study that models the X-ray binary populations from the first 
galaxies of the universe until today. We use as input to our modeling the Millennium II Cosmological Sim- 
ulation and the updated semi-analytic galaxy catalog bv lGuo et al.l (1201 11) to self-consistently account for the 
star formation history and metallicity evolution of the universe. Our modeling, which is constrained by the 
observed X-ray properties of local galaxies, gives predictions about the global scaling of emission from X-ray 
binary populations with properties such as star-formation rate and stellar mass, and the evolution of these rela- 
tions with redshift. Our simulations show that the X-ray luminosity density (X-ray luminosity per unit volume) 
from X-ray binaries in our Universe today is dominated by low -mass X-ray binaries, and it is only at z > 2.5 
that high-mass X-ray binaries become dominant. We also find that there is a delay of ~ 1 . 1 Gyr between the 
peak of X-ray emissivity from low-mass Xray binaries (at z ~ 2.1) and the peak of star-formation rate density 
(at z ~ 3.1). The peak of the X-ray luminosity from high-mass X-ray binaries (at z ~ 3.9), happens ~ 0.8Gyr 
before the peak of the star-formation rate density, which is due to the metallicity evolution of the Universe. 
Subject headings: stars: binaries: close, stars: evolution, X-rays: binaries, galaxies, diffuse background, galax- 
ies: stellar content 



1. INTRODUCTION 

X-ray observations of galaxies are a key component to un- 
raveling the physics of compact object formation. Stellar rem- 
nants and the distribution of the hot interstellar medium are 
best studied via their X-ray emission. With Chandra sev- 
eral global galaxy X-ray correlations hav e been established 
and appear to hold to at least z 
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201 U iLehmer et al.l 12012b . implying a large-scale regular- 



ity to the production of X-ray binaries (XRBs) and hot gas, 
the primary sources of high energy emission in "normal" 
(non-AGN) galaxies. These include a very strong corre- 
lation between total X-ray emission from high-mass XRBs 
(HMX Bs) and the galaxy-wide star-formation rate (SFR) (e.g. 
Ranall iet alj200 3: Gilfa nov et alj2004t|Horn schemei er et al.l 
20051: ILehmer et all 120 lOt iMineoet all 1201 2al) . as well as 
a scaling between total X-ray emission from low-mass 
XRBs (LMXBs) and stellar mass (M„) dGilfanovl l2004t 
Lehm er et al.ll2010t iBoroson et al1l20lTt IZhang et al.ll2012l 
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These correlations present exciting challe nges for models 
of accreting binary evo lution (e.g. Fr agos et al.l 120081 120091: 
iBelczv nski et al] 120081) that predict significant variability in 
these relations based on, e.g., star formation history and 
metallicity. 

Rece nt ultradeep C handra and multiwavelength surveys 
(e.g. [Xue et all 1201 ll) have permitted robust observational 
tests of these correlations in dis tant galaxies (see also the re- 
view of Brandt & Hasinger 2005). In a broad sense, curiously, 
the X-ray emission per unit SFR seems to evolve little with 
redshift, despite the evolutio n in gas metallicity, fo r example, 
over th e sa me epochs (e.g. ILehmer etal.1 120081). iPtak et af] 
(120071) and iTzanavaris & Georgantopoulos (2008) used rich 
multi-wavelength datasets to construct X-ray luminosity func- 
tions (XLFs) of normal galaxies in different redshift bins, and 
found that the XLF shows a statistically significant evolution 
with redshift and that it is the late type galaxies that are driv- 
ing this evolution. 

Despite the significant investment of observing time in X- 
ray studies of distant normal galaxies, and the recent theo- 
retical advances in modeling the X-ray source populations in 
nearby galaxies, our understanding of the cosmological evo- 
lution of populations of compact objects is still in its infancy. 
The first steps in this direction have been made by means 
of semi-analytic, emp irical models (IWhite & Ghoshl 119981: 
Ghosh & White 120011) . Considering a time-dependent SFR, 
these models predicted that the time required for binaries to 
reach the X-ray phase (due to the donor's nuclear evolution 
and angular momentum losses) leads to a significant time de- 
lay between a star-formation episode and the production of 
X-ray emission from X-ray binaries from this populations. 

State-of-the-ar t binary population models (e.g. 

IBelczv nski et al] 120081) provide us with a more physical 
picture of XRB populations as a function of galaxy prop- 
erties. This type of detailed modeling has been applied to 
both starburst and elliptical nearby galaxies, where large 
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populations of individual XRBs are resolved, and detailed 
ages, star-format i on his tories, and metallicities are measured. 
IBelczvnski et all (120041) constructed the first synthetic XRB 
populations for direct comparison with the Chandra observed 
XLFs of NGC 1569, a st ar forming dwarf irregular galaxy. 
iLinden et all (120091 [2010) developed models for the HMXB 
and Be XRBs of SMC, studying the XLF and the spatial 
distribution of this population, and investigating the effect of 
electron-cap t ure supernov ae of massive ONeMg stellar cores. 
Fragos et al. (2008, 2009) performed extensive population 
synthesis (PS) simulations, modeling the two old elliptical 
galaxies NGC 3379 and NGC 4278, constraining the relative 
contribution to the observed XLF from sub-populations of 
LMXBs with different donor and accretor types, and the 
effects of the tra nsient behavior of LMXBs. More recently, 
IZuo & LH (1201 ll) studied the evolution of XRB populations 
in late type galaxies, usin g a modified versio n of the BSE 
population synthesis code (iHurley et al.1 120021) . Considering 
different star-formation history scenarios, they were able 
to derive the X-ray emissivity as well as XLFs of the XRB 
populations. However, their model predictions are not in 
good agreement with observations for most of the Universe 
lifetime: from 4.5 to 13.7Gyr (z= 0-1.4). 

The next step in this long-term effort is to apply these 
models in a cosmological context, in order to understand the 
nature of the accreting populations in high-redshift normal 
galaxies. In this paper, we study the global scaling relation of 
emission from XRB populations with properties such as SFR 
and stellar mass, and the evolution of these relations with red- 
shift. More specifically, we developed a large grid of XRB 
PS models for which the information about the star-formation 
history and metallicity evolution were derived from cosmo- 
logical simulations. Our models, which we have already con- 
strained with observations of XRB populations of the local 
universe, give predictions for the evolution of the aforemen- 
tioned global scaling relations back to the formation of the 
very first galaxies at redshift z > 15. 

The plan of the paper is as follows. In Section 2 we de- 
scribe our simulation tools, the StarTrack PS code and 
Millennium II simulation, and the methodology we follow in 
developing models for the evolution of XRBs on cosmologi- 
cal timescales. Section 3 discusses the necessary bolometric 
corrections we employ in order to directly compare our model 
results with observations. In Section 4, we present the obser- 
vational data we use to constrain our models and the statistical 
analysis we follow. We present the predictions of our maxi- 
mum likelihood models for the evolution of XRB population 
in the high redshift universe in Section 5. Finally, in Section 6 
we summarize the main findings and conclusions of our work. 

2. POPULATION SYNTHESIS MODELING 

2.1. Population Synthesis Code: StarTrack 

The main t ool we use to perform the PS simulations is 
StarTrack (IBelczvnski et al J 120021 12008). a state-of-the- 
art code with special emphasis on processes leading to the for- 
mation and further evolution of compact objects. Both single 
and binary star populations are considered. The code incorpo- 
rates calculations of all mass-transfer phases, a full implemen- 
tation of orbital evolution due to tides, as well as estimates of 
magnetic braking. StarTrack has been extensively tested 
and calibrated using detailed mass-transfering binary star cal- 
culations and observations of binary populations. For a com- 
prehensive summary of StarTrack and its model parame- 



ters see IBelczvnski etai](l2008l) 

Recently the StarTrack code has undergone three ma- 
jor revisions. The first one allows for the updated stel- 
lar winds and their re-c alibrated dependence on metallicity 
(IBelczvnski et al.l[2010l) . This revision is fully incorporated 
into our results. The second update, based on the fully con- 
sistent supernova simulations, resulte d in the revised neutron 
star and black h ole mass spectrum (IBelczvnski et al.l 120121; 
iFrver et aT]|2012l) . Finally, the most recent upgrade involves a 
more physical treatment of donor stars in common envelope 
events via usage of actual value of A parameter for whic h usu- 
ally a constant value is assumed (Dominik et al. 2012). The 
last two code revisions are not employed in our simulations as 
these were performed well before the revisions became avail- 
able. 

Several XRB P S stu die s that used Sta rT rack 
(IBelczvnski et all 12004 IBelczvnski & Taarnl 12004 
Fragos_etaL| |2008, 2009T 120101: Linden et all " 120091 12010) 
have been pivotal in interpreting observations of nearby 
galaxies and showed that despite the significant number of 
free parameters (~ 10) and the complexity of the physical 
processes involved, PS models can provide robust results 
that can be used to understand the properties of more distant 
galaxies. 

We should stress here that our PS code only models XRBs 
formed via the evolution of primordial isolated binaries, i.e. 
the field XRB population. A dynamically formed population 
of LMXBs can have a significant contribution to the integrated 
X-ray luminosity of some globular cluster rich elliptical 
galaxies, where in certain cases more th an half of the bright 
LMXBs reside in globular clusters (e.g. Hump hrey & Buotel 
2008). The detailed modeling of the dynamically formed 
LMXBs, although possibly important for this subset of ellip- 
tical galaxies, is outside the scope of this paper, and in fact it 
is not feasible given the tools that exist to date for simulations 
of such scale. The only PS studies to date of dynamically 
formed LMXBs have been done for a few specific Galactic 
globular clusters and only in order to identify t he various for- 
mation channels a nd their relative importance (llvanova et alj 
2006, 2008L l2010h . 

In the work presented here we consider the global XRB 
population, i.e. XRBs in a general galaxy population of both 
early and late galaxies, where the contribution of dynamically 
formed LMXBs is much smaller. In the low redshift Uni- 
verse appr oximately half of the stellar mass is in early type 
galaxies dBell et al.l l2003). Furthermore, approximately half 
of early type galaxies are old, globular-cluster rich ellipti- 
cals, and finally in these globular-cluster rich ellipticals ap- 
proximately half of the bright LMXBs are located in globu- 
lar clusters. This effectively translates to a possible error of 
10% -25% in the determination of the integrated X-ray lu- 
minosity coming from the LMXB population when neglect- 
ing the dynamically formed population. This contribution of 
the dynamically formed population is decreasing further as 
we move to higher redshifts. However, comparing the X- 
ray luminosity per un it stellar mass com ing from LMXBs in 
star-forming galaxies (iLehmer et al.l2010l) . where the dynam- 
ically formed LMX B population is negligible, and in ellipti- 
cal galaxies dBoroson et alJl2011l) . shows that any correction 
that one should make due to the contribution of dynamically 
formed LMXBs is smaller than the observational uncertain- 
ties. 

2.2. The Millennium-II Simulation 
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The Millennium-II Simulation (Bovl arTkolchin et alj 
2009) is a very large N-body simulation of dark matter 
evolution in the concordance ACDM cosmology. The sim- 
ulation uses ~ 10 10 particles in a lOQMpc/h box, resulting 
in a spatial resolution of lkpc/h and a mass resolution 
of 6.89 x 1O 6 M //i. 67 full snapshots of the simulation 
were stored, with time intervals between them varying from 
~ 30Myr to ~ 300 Myr, which were then used to derive 
merger trees of the dark-matter halos. 

Using the stored h alo/subhal o merg er trees of the Millen- 
nium II simulation, iGuo et all d201 lh recently r e leased an 
updated semi-analytic galaxy catalog. IGuo et af] d201 lh re- 
vised the modeling of several physical processes involved in 
galaxy formation, such as gas stripping and supernova feed- 
back. This newly released catalog yields an excellent fit to the 
galaxy mass and luminosity function of low-redshift galaxies 
over five orders of magnit ude in stel l ar ma ss and nine mag- 
nitudes in luminosity. The lGuo et alj d201 lb catalog provides 
a wealth of information for each galaxy in the cosmological 
simulation, including galaxy type, SFR, metallicity, hot and 
cold gas content, central black hole (BH) mass, and rest-frame 
absolute magnitudes in the SDSS filters, as a function of time. 

As already mentioned, the Millennium II simulation is a 
dark matter only simulation, and the baryonic matter is mod- 
eled as an add-on, using semi-analytic prescriptions. De- 
spi te the high comp lexity of the galaxy formation modeling 
by IGuo et all (1201 ll) . which takes into account most physi- 
cal processes involved, it remains an approximate model. In 
that sense, the detailed evolution of a specific galaxy formed 
in the simulation may not be fully con sistent with more 
advan ced galaxy formation simulations (IScannapieco et all 
20121). Howeve r, the semi-analytic prescriptions used by 



Guo et all (1201 lh have been calibrated so that the statistical 
properties (e.g. mass-function, global star-formation history, 
stellar metallicity distribution) as a function of redshift of 
the modeled galaxy population are consistent with available 
observational data. Figure Q] sh ows the p r edict ions of the 
Millennium-II simulation and the lGuo et alj d201 lh modeling 
for the SFR density (top panel), the stellar mass density (mid- 
dle panel), and the metallicity of the newly formed stars as 
a function of redshift. For comparison observed data for the 
SFR density and stellar mass density are also plotted (upper 
and middle panel respectively), showing a good agreements 
between model and observations. 

2.3. The Millennium-II simulation as Initial Conditions of 
Populations Synthesis Models 

The first step of our analysis is to derive for each combi- 
nation of PS model parameters, and each value of metallicity, 
the evolution of the specific X-ray luminosity (X-ray lumi- 
nosity per unit stellar mass) of a coeval population of XRBs 
as a function of its age. Assuming a single instantaneous star- 
burst, we follow the evolution of the XRB population, keeping 
track of the total (cumulative) X-ray luminosity of the whole 
XRB population, and the contribution of the different sub- 
populations (e.g. XRBs with different types of accretors) as 
a function of the age of the population. We note that in this 
work we define as a LMXB, a XRB with donor star whose 
current mass is less than 3M Q , while a HMXB is any XRB 
with a donor star more massive than 3 M Q . 

In the top panel of Figure [2] we show the specific X-ray 
luminosity of XRBs, X-ray luminosity per unit stellar mass 
(Lx/M*), versus age for a single coeval population of bina- 
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Figure 1. The top panel show s the SFR density evolution with redshift, 
as predicted by the Guo et al. 120111) modeling (solid line). Observa- 
tional determinations of the SF R density at different redshifts (grey squares) 
are shown for compari son (Schiminovich et al. 2005; Bouwens et alj |2007t 
IReddv & Steidel 2009). The stellar mas s density as a function of redshift, 
as predicted by the the [Guo et al. ( 2011) modeling, is shown in the middle 
panel (solid line) , as we ll as a compilation of observed data (grey squares) by 
Marchesini et al. (2009). The bottom panel shows the mass weighted mean 
metallicity of the newly formed stellar mass as a function of redshift. The 
dark grey area shows the metallicity range with which 68.2% of the new stars 
form, as a function of redshift, while the light grey area corresponds to the 
95.4% of the newly formed stars. 

ries at solar metallicity. These curves were generated using 
our "reference" model (Model 245; see also sections 2.4 and 
4.2). In the bottom panel of Figure [2] we show how these 
curves vary if instead initial metallicities of ~ 0.1Z Q (green 
dashed) and ~ 1.5Z Q (magenta dot-dashed) are chosen. As 
seen in the top panel of Figure [2] over the lifetime of the stars 
produced in a given star-formation event, there is an initial 
ramp-up in the formation of HMXBs over the first ~ 5 Myr as 
the first compact objects form, followed by a period of HMXB 
dominance that lasts for ~ 100- 300 Myr. At this point, the 
LMXB population emerges, peaks, and then slowly declines 
by ~ 1-2 orders of magnitude out to ~ 14Gyr timescales. In 
the lower panel of Figure [2] we show the X-ray binary lumi- 
nosity evolution with respect to solar for low and high metal- 
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Figure 2. The evolution of an XRB population formed from a single-burst 
population of stars. Top panel: Specific bolometric X-ray luminosity, i.e. lu- 
minosity per unit stellar mass, of a coeval stellar population as a function of 
the population's age for our "reference" model (Model 245; see also sections 
2.4 and 4.2) and for solar metallicity, broken down into HMXB (blue dotted) 
and LMXB (red long-dashed) contributions. After an initial ~ 5Myr ramp- 
up, HMXBs dominate the X-ray power output for ~ 100-300Myr. Fol- 
lowing this, LMXBs take over, peak in formation at ~ 0.5 - 1 Gyr, and then 
passively fade with increasing age. The bumps and wiggles present at > 4 Gyr 
are due to statistical fluctuations as bright sources turn on and off. Bottom 
Panel: X-ray binary evolution for the cases of metallicities Z = 0. IZq (green 
short-dashed) and Z = 1.5Z0 (magenta dot-dashed) compared to solar (i.e., 
L x [Z]/Lx[Z @ ]). 

licity cases. We note that the X-ray luminosity is expected to 
be significantly higher for low-metallicity regions, due to the 
fact that the young massive O/B stars, which are the precur- 
sors to compact object accretors in X-ray binaries, will lose 
less stellar mass from line-driven winds over their lifetimes 
for lower metallicities. This results in more numerous and 
more massive BH populations, and subsequently more lumi- 
nous X-ray binary populations (see also section 5.1 for a de- 
tailed discussion). Some evidence has already indicated that 
the X-ray luminosity from XRBs per unit of SFR (Lx/SFR) 
for low-metallicity dwarf gal axies may be enha nced com- 
pared to more typical galaxies (Kaaret et al. 201 1). 

The Millennium II simulation and the iGuo et"afl (1201 lh 
semi-analytic galaxy catalog provide us information about the 
evolution of the properties of each galaxy in the simulation. 
In this first study, we consider only the global characteris- 
tics and scaling relations of the XRB evolution across cosmic 
time. Having this in mind, we consider the whole simulation 
box of the Millennium II simulation as representative of the 
Universe, and we keep track of the total new stellar mass that 
is formed between each snapshot of the simulation. In addi- 
tion, we keep track of the metallicity of the new stellar mass 
formed. In practice, we consider 9 metallicity bins, centered 
at the values reported in Table [1] and calculate how much new 
stellar mass was formed at each snapshot and at each metallic- 
ity bin. Finally, we assume that the new stellar mass is formed 
with a constant SFR between each snapshot. This results in a 
self-consistent prescription of the star-formation history and 
metallicity evolution of the Universe. 



We now have all the necessary information in hand to con- 
struct models for the evolution of XRBs across cosmic time. 
We convolve the relations we derived in the first step, with the 
star-formation history for different metallicity ranges. The re- 
sult is a mixed XRB population, with a broad range of ages 
and metallicities, representative of the XRB population of the 
Universe, at each redshift. From this convolution we derive 
model predictions such as the X-ray luminosity from XRBs 
per unit volume, stellar mass, or SFR, as a function of red- 
shift. 

2.4. Parameter study 

StarTrack PS models have a significant number of free 
parameters, which can be categorized into two groups. In the 
first group there are parameters that usually correspond to the 
distributions of the initial properties of the binary population, 
such as the initial mass function (IMF), the initial binary mass 
ratio (q), and the distribution of initial orbital separations. For 
these parameters, we primarily use fixed values derived from 
the most updated observational surveys. Only in the case that 
there is an ongoing debate regarding a parameter we actually 
treat it as a free parameter. In the same group one should add 
parameters such as the stellar wind prescriptions, and the natal 
kick distribution of neutron stars (NSs), as these parameters 
are also derived from observational data. The second group 
is that of the few "true" free parameters that correspond to 
poorly understood physical processes which we are not able 
to model in detail. In this second group belong parameters 
such as those involved in the common envelope (CE) phase 
modeling. 

We created the largest ever grid of 288 PS models. Each 
of these 288 models was run at 9 different values of metal- 
licity (Z= 10 _4 -0.03), where we followed the evolution of 
~ 45 million binaries per model. The whole set of simula- 
tions required ~ 2 million CPU hours of computational time. 
In this grid, we varied 6 parameters (i n addition to the metal- 
licity) known from previous studies (IBelczvnski et "aLI I2007L 
120 lOt iFragos et al.ll2008l 120 lOt iLinden et al.1 120091) to affect 
the evolution of XRBs and the formation of compact objects 
in general. Namely, we varied the CE efficiency (o.ce x A = 
0. 1-0.5), the IMF (Kroupa 2001 or Kroupa & Weidner 2003), 
the initial mass ratio (q) distribution (flat with q = 0- 1, twin 
with q = 0.9- 1 .0, or a mixture of 50% flat and 50% twin), the 
stellar wind strength (rfaM = 0.25, 1 .0, or 2.0; parameter with 
which we multiply the stellar wind prescription described in 
iBelczvnski et al.ll2010l) . and the distribution of natal kicks for 
BHs formed though direct core collapse (zero kicks, or 10% 
of the standard kicks for NS). Finally, we allow for various 
outcomes of CE phase, in particular taking into account the 
possible CE inspirals with Hertzspmng gap donors that termi- 
nate binary evolution b arring the subsequent XRB formation 
(Belczvnsk i et al.ll2007l) . The grid of 288 models we created 
includes all the possible combinations of the parameters we 
varied. The various parameter values used are shown in Ta- 
bleQ] 

We note that in our calculations we combine o:ce and A 
into one CE parameter, where A is a measure of the central 
concentration of the donor and the envelope binding energy. 
In the rest of the text, whenever we mention the CE effi- 
ciency cecE, we refer i n practice to the product cuce x A (see 
IBelczvnski et al.ll2008l for details). 

3. BOLOMETRIC CORRECTIONS TO THE CHANDRA X-RAY 
BANDS 
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Table 1 

Model Parameters 



Parameter 



Notation 



Value 



Reference 



Initial Orbital Period distribution F(P) 

Initial Eccentricity Distribution F(e) 

Binary Fraction /j,j n 
Magnetic Braking 

Metallicity Z 
IMF (slope) 

Initial Mass Ratio Distribution F(q) 

CE Efficiency oce 

Stellar wind strength Vv/ind 

CE during HG 

SN kick for ECS/AIC NS 

SN kick for direct collapse BH 



flat in logP 
Thermal F(e) ^ 
50% 



0.0001,0.0002, 0.005,0.001, 
0.002, 0.005,0.01,0.02, 0.03 
-2.35 or -2.7 
Flat, twin, or 50% flat + 50% twin 
0.1,0.2, 0.3, or 0.5 
0.25, 1.0, or 2.0 
Yes or No 
20% of normal NS kicks 
Yes or No 



Abt(1983) 
Heggie (1975) 

Ivanova & Taam (2003) 



Kroupa (2001); Kroupa & Weidner (2003) 
Kobulnickv & Fryer (2007); Pinsonneault & Stanek (2006) 
Pods iadlowski e t al. (20031 
Belczvnski et al. (2010) 
Belcz vnski et al. (2007) 
Linden et al. (2009) 
Fragos et al. (2010) 



a Electron Capture Supernova / Accretion Induced Collapse 

StarTrack, our PS synthesis code, keeps track of the 
mass-transfer rate as a function of time for every mod- 
eled XRB. From this mass-transfer rate, we derive the bolo- 
metric luminosity based on the prescriptions presented by 
iFragos et all (120081 12009b . which also take into account the 
transient behavior of XRBs. However, in order to compare 
our models with observational data, we need to estimate the 
X-ray luminosity in the specific energy band where the data 
w ere taken. 

iMcClintock & Remillardl (120061) and lWu et all feOlO) com- 
piled two samples of RXTE observations of Galactic NS and 
BH XRBs at different spectral states, for which they also cal- 
culated the best fit parameters of simple spectral models. Out 
of the two lists (which in some cases contained multiple ob- 
servations of the same object at different spectral states) we 
chose all XRBs for which there was a distance measurement 
and a BH mass determination for BH XRBs. For the case 
of NS XRBs, we assumed a nominal mass of 1.4M©. Us- 
ing the appropriate XSPEC models, we calculated bolometric 
correction factors between the bolometric X-ray luminosity 
(0.03- 100 keV) inferred by the best fitting model for each ob- 
servation, and the X-ray luminosity in a specific energy band. 
We should note here that the fraction of the energy emitted 
by XRBs outside the 0.03- 100 keV range is negligible, and 
hence the terms bolometric luminosity and X-ray bolometric 
luminosity can be used interchangeably. For each system we 
calculate these bolometric correction factors with and with- 
out absorption corrections. These factors would correspond 
to the observed and intrinsic X-ray luminosity at a specific 
energy band respect ively. 

According t o IMcClintock & R emillard (2006) and 
IWu et al.l (120101) the XRBs in their compilations are in either 
a low-hard spectral state, where the spectrum is dominated 
by a power-law component, or in a high-soft state, where the 
spectrum is dominated by the t hermal emission of the disk. 
IMcClintock & Remillardl (120061) also use a third category of 
a "steep power-law" state for BH XRBs with accretion rates 
close to the Eddington limit. These systems do show a soft 
spectrum, which however can be fit better with an absorbed 
steep power-law than with a thermal disk model. For our 
purposes, we do not treat separately systems in the high-soft 
and steep power law states, as they turn out to have very 
similar bolometric correction factors. 

In our calculations we treat separately NS and BH 
XRBs. The transition between low-hard and high-soft 
state in both cases happens at a luminosity of ~ 5%LEdd 
(IMcClintock & R emillard 2006). Figure [3] shows the correc- 
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Figure 3. Bolometric corre ction factors for the 0.3 -8keV band, for each 
of the observations from the McClintock & Remillard (2006) and Wu et al. 
1 2010) samples. Star symbols denote systems in the low-hard state, while 
open squares systems in the high-soft state. Blue color corresponds to NS 
XRBs and red to BH XRBs. The dashed vertical line shows the approximate 
Eddington ratio at which the state transition happens. The horizontal lines 
denote the mean bolometric correction of NS or BH XRBs at the low-hard or 
high-soft state respectively. 



tion factor (j] Bo i) between the bolometric X-ray luminosity 
and the intrinsic (after correcting for absorption) luminosity 
in the 0.3-8 ke V band (i.e. r\Boi^Boi = £o.3-8kevX as a function 
of the Eddington ratio of the source. For each energy band 
we calculate the mean and the variance of the bolometric cor- 
rection for the high-soft and low-hard states, and BH and NS 
XRBs separately. Furthermore, we perform our calculations 
with and without correcting for absorption due to circumstel- 
lar or interstellar gas, corresponding to the observed and the 
intrinsic luminosity accordingly. The full list of the calculated 
bolometric corrections for different bands can be found in Ta- 
ble [2] The mean values of the calculated bolometric correc- 
tions f or the different ene rgy bands are used in equations (1 )- 
(6) of iFragos et al.l Q2008) in order to compare our PS models 
with observational data. The variances we calculate can en- 
ter in the estimation of the uncertainties of our models. In 
the rest of the paper, we use the absorption uncorrected val- 
ues in our comparison, as these may be most easily compared 
with observed X-ray luminosities of galaxies. This is effec- 
tively equivalent with the assumption that all galaxies in the 
universe have similar intrinsic absorption as the Milky Way. 

Although these bolometric corrections are based on spec- 
tral fits of RXTE data which have a low-energy limit of 
~ 2.5 keV, the RXTE band covers most of the Wien tail 
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of the disk black-body component, allowing a reliable mea- 
surement of the disk temperatu re (e.g. Figures 4.8, 4.9, in 
iMcClintock & R emillard 200<|. A comparison between the 
disk temperature based on RXTE measurements, with sim- 
ilar measurements for the same objects in similar states de- 
rived from spectra extending to lower energies (using data ob- 
tain ed with XMM-Newton, Suzaku, Beppo - SAX, or ASCA; 
e.g. iKotani et alj 120001: iKubota et alj|200l: iGou et alj 120091: 
Caball ero-Garcia et a l. 2009; Tamu ra et al. H2012T shows very 
good agreement, indicating that indeed the relatively re- 
stricted soft X-ray coverage of RXTE does not significantly 
bias the temperature of the thermal component of sources in 
the high-state. 

4. OBSERVATIONAL CONSTRAINTS AND STATISTICAL MODEL 
COMPARISON 

4.1. Available observational constraints 

To constrain our model parameters, we compare our mod- 
els to observations of the X-ray properties of local galax- 
ies. There have been several recent studies, where high 
quality multi-wavelength data are used to derive scaling re- 
lations between total X-ray emission from HMXBs and the 
galaxy-wide SFR as well as a scaling between total X-ray 
emissio n from LMXBs and stellar mass, in the local uni- 
verse. Boroso n et al.l (1201 ll) selected a sample of 30 nor- 
mal early-type galaxies, for all of which optical spectroscopy 
and Chandra data were available. From the optical data 
they estimated the stellar mass, while from spectral fitting 
of the X-ray data they estimated the contribution of unre- 
solved LMXBs to the diffuse X-ray emission . We aver- 
aged the measurements for the Boros orTet al.l (120111) sam- 
ple to calculate the X-ray emission from LMXBs per unit 
stellar mass, fin ding log(L x ,lmxb s [0.5 -8 keV] /M« ) = 29.18 ± 
0.17ergs" 1 MQ 1 . ILehmer et all (120101) used Chandra obser- 
vations of star-forming galaxies to study the galaxy-wide 
XRB emission and its correlation with SFR and stellar mass. 
They assumed that X-ray emission in the hard X-ray band 
(2- lOkeV) of these galaxies can be decomposed into two 
components. The first component is related to HMXBs and 
should be proportional to the SFR of the galaxy, while the 
second component is related to LMXBs and should be pro- 
portional to the stellar mass of the galaxy. This simple model, 
which is described by the equation L x [2-iokeV] = oM^ + (3SFR, 
was fitted to their data, and the best fitting values were 
found to be: loga = log(L x , L MXB S [2-mk e v]/M f ) = 28.96 ± 
0.34 erg s _1 Mq' and logfi = log(L x , HMXBs[2 -w ktV]/SFR) = 
29.21 ± 0.34 erg s^M^ 1 yr. Most recently, iMineo et all 
d2012al) . based on a homogeneous set of X-ray, infrared and 
ultraviolet observations from the Chandra, Spitzer, GALEX 
and 2MASS archives, studied populations of HMXBs in a 
sample of 29 nearby star-forming galaxies and their relation 
to the SFR. In agreement with previous results, they found 
that HMXBs are a good tracer of the recent star-formation 
activity in the host galaxy and their collective luminosity 
scales with the SFR as log(L x ,HMXB S [0.3-8keV]/SFR) = 29.47 ± 
0. 4ergs~'M^ 1 yr. 

ITzanavaris & Geor gantopoulosl (120081) used multi- 
wavelength datasets of the CDFS, E-CDFS, CDFN, and 
XBootes fields to construct XLFs of normal galaxies in 
different redshift bins. The integral of these XLFs gives us 
the total X-ray lumino sity per unit volume of normal galaxies 
in three redshift bins. ITzanavaris & Geo rgantopo ulosl (120081) 
used the soft, 0.5-2keV, Chandra energy band, as Chandra 



is most sensitive in this band. However, this band likely 
contains also significant emission from hot diffuse gas. 
Hence we can only use these observed values as upper 
limits when comparing them to our simulations which do 
not include any emission from hot gas. We should note 
here, that in the last years there have been several X-ray 
stacking stud i es, es p ecially in the Chandra deep fie lds (e.g. 
ILehmer et al.1 120071 12001 ISvmeonidis eTail l20Tlh . where 
the X-ray luminosity per unit stellar mass or unit SFR, for 
early and late type galaxies respectively, has been measured 
at z » 0— 1.4. However, a proper comparison of our models 
with these observational surveys would require a careful 
matching of the sample selection of these studies, which is 
outside the scope of this paper. 

4.2. Model likelihood calculation 



We compa r e the results by iBoroson et al.l (12011b and 
Lehm eret al.l (|2010) to the prediction of our models for 
Lx.r.MXB.jM* at z = and in two diff erent energy bands, a nd 
the results by lMineo et al.l (12012a!) and lLehmer et al.l (120101) to 
our models' predicted L x ^mxbs/SFR, again in two different 
energy bands. Furthermore, we compare our model predic- 
tions for th e X-ray luminosity per unit vo l ume w ith values 
reported by Tzanavaris & Georgantopoulos (2008), which we 
use as upper limits. The quantitative comparison is done by 
calculating the likelihood of the observations given a model, 
for each model, and then comparing the likelihood values. 

Assuming that the reported errors for all observed quantities 
are gaussian in the log space, then the likelihood of the data 
given a model becomes: 



Li(D\M,fj)= \[ GiMijJjfXj^Pifj) 



(1) 



7=1,7 



where G is the Gaussian function, fij and aj for j = 1,4 
are the observed values and the associa ted errors at 
z = for los(L Y jMYR.m s-kv. WMJ (IBoroson et all 
1201 ll) . los(L Y fMYBrn-mi.Pvi /MJ (ILehmer et alj 120101). 
log{L XflMXBs [ 2 -\Qk e v]/SFR) (ILehmer et all 120 id) , and 
log(L x , HMXBs[0 3-&k e v]/SFR) (Mine o et al]l2012al) . /x, and aj 
for j = 5,7 are the observed values and the associated errors at 
z ~ 0.13, z ~ 0.38, and z ~ 0.78 for log (L x [o.5-2keV] /Volume) 
(ITzanavaris & G eorgantopoulos 2008f. M,j are the predic- 
tions of the i th model for each of these quantities and fj is the 
fraction of the luminosity density /!_,■ assumed to be due to 
X-ray binaries. P(fj) is the prior probability for a given value 
of fj. Since we are not attempting to constrain the X-ray 
binary fraction / we marginalize over /, i.e., 

LijiMij \D) = J I.Al) \M U , fj)P( f j)dfj (2) 

For the z = points we assume that the luminosity densities 
are purely due to X-ray binares and therefore /=1. In this 
case P(f j = \,4) = S(fj =1). We assume a flat prior for the X-ray 
binary fraction for the z > luminosity densities, i.e., P(fj) = 
Cj, where the constants Cy are chosen so that the likelihood 
function remains normalized to unity. Noting that: 



G(M L j-fj Nl aj)Cjdfj = [" G(Mij- n' j^Cjdfi'j = 
Jo 



(3) 
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Table 2 

Mean value and standard deviation of the bolometric correction (r?j, D ;) factor at different energy bands. 



Energy band 



Absorption Corrected 
High-Soft State Low-Hard State 

NS BH NS BH 



Absorption Uncorrected 
High-Soft State Low-Hard State 

NS BH NS BH 



0.3-2keV 

0.3-7keV 

0.3-8keV 

0.3-10keV 

0.5-2keV 

0.5-10keV 

2-10keV 



0.06±0.09 
0.61 ±0.09 
0.67±0.09 
0.77±0.08 
0.06±0.09 
0.77±0.08 
0.71 ±0.08 



0.28±0.20 
0.75±0.18 
0.78±0.16 
0.82±0.13 
0.26±0.19 
0.80±0.11 
0.54±0.12 



0.08±0.06 
0.29±0.12 
0.31±0.12 
0.36±0.13 
0.08±0.06 
0.36±0.13 
0.28±0.09 



0.13±0.11 
0.36±0.13 
0.39±0.13 
0.43±0.13 
0.12±0.09 
0.42±0.13 
0.30±0.11 



0.05±0.08 
0.57±0.09 
0.63±0.09 
0.73±0.08 
0.05±0.08 
0.73±0.07 
0.68±0.07 



0.09±0.10 
0.48±0.14 
0.51±0.13 
0.55±0.12 
0.09±0.10 
0.55±0.12 
0.46±0.09 



0.05±0.03 
0.25±0.08 
0.27±0.09 
0.32±0.10 
O.O5±0.O3 
0.32±0.10 
0.27±0.08 



0.08±0.06 
0.30±0.08 
0.33±0.10 
0.38±0.10 
0.07±0.05 
0.37±0.10 
0.30±0.11 



Table 3 

List of selected PS models that appear in the text and figures of this paper. A complete list of all models can be found in the 

online version of the article. 



Model 




IMF exponent b 


Vwind C 


CE-HG d 


q distribution 6 




Ranks 


log(Likelihood ratio) 11 


245 


0.1 


-2.7 


1.0 


No 


50-50 


0.1 


1 


0.0000000 


229 


0.1 


-2.7 


2.0 


Yes 


50-50 


0.0 


2 


-0.035048515 


269 


0.1 


-2.7 


1.0 


Yes 


50-50 


0.1 


3 


-0.11226917 


205 


0.1 


-2.7 


2.0 


No 


50-50 


0.0 


4 


-0.13235138 


249 


0.1 


-2.35 


2.0 


No 


50-50 


0.1 


5 


-0.28666705 


273 


0.1 


-2.35 


2.0 


Yes 


50-50 


0.1 


6 


-0.28926393 


53 


0.1 


-2.7 


1.0 


No 


flat 


0.1 


20 


-1.3945094 


246 


0.2 


-2.7 


1.0 


No 


50-50 


0.1 


25 


-2.0332112 


241 


0.1 


-2.35 


1.0 


No 


50-50 


0.1 


31 


-2.3582808 


253 


0.1 


-2.7 


2.0 


No 


50-50 


0.1 


14 


-0.90617906 


261 


0.1 


-2.7 


0.25 


No 


50-50 


0.1 


43 


-3.4548391 


197 


0.1 


-2.7 


1.0 


No 


50-50 


0.0 


13 


-0.80759106 



a CE efficiency parameter 

b Exponent of the high-mass power law component of the IMF:|Kroupa 1 2001 !) (-2.3 5) orlKroupa & Weidnei 1 2053) (-2.7). 
c Stellar wind strength parameter with which the "standard" ( Belczynski et al. 2010) stellar wind prescription is multiplied. 



Yes: all possible outcomes of a CE event with a HG donor are allowed, No: A CE with a HG donor star will always result to a merger. 
e Binary mass ratio distribution. "50-50" indicates half of the binaries originate from a "twin binary" distribution and half from flat mass ratio 
distribution. 

^ Parameter with which the "standard" Hobbs et al. 12005) kick distribution is multiplied for BHs formed though a SN explosion with negligible 
ejected mass. 

° Ranking of the model in our statistical comparison. 

b Ratio of the likelihood of the observations given a model to our maximum likelihood model (log L(0\Mi) - \og HO\M re f)). 
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Figure 4. Ratio of the likelihood of the observations given a model to our maximum likelihood model. 



8 



T. Fragos et al. 



The marginalized likelihood is then 
I-i (Mi n> = 



(at 



(4) 



•Mi. 



M 



2a) 



= ~[[G(Mij;nj,aj)x~[[ erf KI ~\"' '' J ' -erf 
7=1,4 ;=5,7 

Figure|4]shows the likelihood ratio of each of our 288 mod- 
els to the maximum likelihood model (L(D\M re f J). This ef- 
fectively means that we renormalized the likelihood of the 
observations given a model, so that our maximum likelihood 
model has a likelihood value of 1. A list of the likelihood 
values of each model can be found in Table 3] One should 
note that the scale of the Y-axis in Figure |4] is logarithmic, 
and only a small fraction of all the 288 models we consid- 
ered have a likelihood close to the our reference, maximum 
likelihood model (Models 245). Namely, only 14 out of 288 
models have likelihood values within a factor of 10 from the 
maximum likelihood model, and only 24 within a factor of 
100 respectively, which is an encouraging evidence that our 
calculated likelihood is a good discriminator between the dif- 
ferent models. One other apparent feature of Figure 2]is that 
models Models 100-199 have systematically lower likelihood 
values than the rest. In all these models the initial binary mass 
ratio is assumed to be "twin", i.e. a flat distribution between 
0.9 and 1.0, which always results in an LMXB population in- 
consistent with observations (see also Section 5.1 for further 
discussion). 

Figure [5] shows the X-ray luminosity in the 2 - lOkeV 
range, from the LMXB population, per unit stellar mass (top 
left panel) and X-ray luminosity in the same range, from 
the HMXB population, per unit SFR (top right panel) as 
a function of redshift. The six models shown in this fig- 
ure are the ones that satisfy within one a all the observa- 
tional constraints, which are also the six models with the 
highest likelihood values. The observed values of these 
two quantities for z = bv iLehmer etai] (|2010l) . and by 
IBoroson etafl d2011l) and lMineo et al.1 d2012al) after convert- 
ing them to the 2- lOkeV range, are also shown for com- 
parison. The bottom panel of Figure [5] shows the X-ray lu- 
minosity density coming from the the whole XRB popula- 
tion as a f unction of redshift. The observational constraints 
derived by Tzanavaris & Georgantopoulos (2008), shown as 
green squares in the figure, are used as upper limits in our sta- 
tistical analysis. Figure ^demonstrates the robustness of our 
model predictions. All six maximum likelihood models ex- 
hibit a very similar behavior for the whole range of redshifts, 
and the differences between the different models at higher 
redshifts are equivalent to the uncertainties of the observed 
values in the local universe. 

4.3. The specific X-ray luminosity of local elliptical galaxies 

iKim & Fabbianol ((2010), using a sample of local elliptical 
galaxies with measured stellar ages, found that the fraction of 
bright LMXBs in "young" elliptical gal axies is higher com - 
pared to "old" ellipticals. More recently, Zhang e t al.l (120121) . 
using a similar sample of nearby ellipticals, reported a pos- 
itive correlation between the normalization of the XLFs of 
these galaxies, and the luminosity average age of their stellar 
population, as estimated by single stellar population spectral 
energy distribution fitting results. The specific X-ray lumi- 
nosity in these galaxies (Lx/M*), which depends both on the 
shape and the normalization of the XLF, shows little to no evo- 
lution with the estimated stellar age of the galaxies (see also 



IBoroson et al.ll201lT) . which at face value seems to be incon- 
sistent with our prediction for the evolution of Lx/M* with 
redshift. However, the galaxy sample in both of these stud- 
ies shows also a strong correlation of the estimated stellar age 
with the globular cluster specific frequency (Sn), as typical 
nearby old giant ellipticals tend to have a high 5V- At the same 
time a strong correlation between Lx /M* and Sm has been ob- 
served in these galaxies (IBoroson et alj20lTb . In addition, the 
age estimates of the observed galaxies are derived under the 
assumption of a single age stellar population, and assigning 
single ages to "young" ellipticals will suffer from large uncer- 
tainties in the stellar ages due to very youn g populations mix - 
ing with old stellar populations (see, e.g.. Ildiart et al.l f2007). 
Hence, it is hard to disentangle whet her the observed corre- 
lation reported by Zhang et al.l (12012b means that the specific 
X-ray luminosity of the LMXB population in these galaxies 
indeed does not evolve with stellar population age, or it is 
an effect of the underlying correlations in the properties of 
the galaxy sample they used, and uncertainties in the age es- 
timates. In either case, we should stress again here that our 
models consider only field XRBs and a general galaxy popu- 
lation, where the dynamically formed XRBs are expected to 
have minimal contribution. In contrast, the observed correla- 
tions discussed above are driven by globular cluster rich old 
elliptical galaxies whose LMXB population i s probably d om- 
inated by dynamically formed LMXBs (e.g., Irwin 20051). 

Figure [6] shows the evolution of the specific X-ray lumi- 
nosity from a coeval stellar population (similar to Figure [2] 
but in the 0.3-8keV energy range) as a function of the stel- 
lar population age for our "reference" model (Model 245; see 
also sections 2.4 and 4.2) and for solar metallicity. In the 
same figure we have over-plotted the LxjM* in the same en- 
ergy rang e for elliptical galaxi es with stellar age estimates 
from the Boroson e t al.l (1201 ll) sample. The galaxy sam- 
ple in Boroson et al.l (1201 ll) is very similar to the sample by 
Zhan g etal.1 (I2012IK We want to stress here that a direct com- 
parison of the model curve to the observed data points is not 
appropriate for two main reasons: the model curve neglects 
the dynamically formed LMXB population which can be sig- 
nificant for some of these galaxies, and the age estimates of 
the observed galaxies are derived under the assumption of 
a single age stellar population. However, assigning single 
ages to young ellipticals will suffer from large uncertainties 
in the stellar ages due to very "young" populations mixing 
with old stellar populations. Keeping in mind these two im- 
portant caveats, we can make a few instructive observations: 
(i) Comparing the model curve to the observed LxjM* of 
the "old" (~ 9— 12Gyr) ellipticals, we find that our models 
under-predict the LMXB specific luminosity. This is to be ex- 
pected, as in these galaxies the dynamically formed LMXB 
population is significant, and our models consider only the 
field LMXB population, (ii) Comparing the model curve to 
the ellipticals with ages ~ 4-5Gyr (which also tend to have 
lower specific globular cluster frequencies), we find that our 
models predict the same Lx /M* as the observations. This can 
be also seen in the top left panel of Figure [5] The average 
mass-weighted age of the stellar population at z = based on 
the Millennium II simulation is ~ 4.7 Gyr, and the predicted 
from our models Lx /M* (coming from LMX Bs) is in very 
good a greement with the Lx/M* as reported by lBoroson et al.l 
d201 ll) . (Hi) Our model curve is in disagreement with the 
observationally determined Lx /M* of the "young" ellipticals 
(ages of ~ l-2Gyr). However, these are the galaxies for 
which the single age stellar population assumption most likely 
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Figure 5. X-ray luminosity in the 2- lOkeV range, from the LMXB population, per unit stellar mass (top left panel) and X-ray luminosity in the same 
range, from the HMXB population, per unit SFR (top right panel), as a function of redshift, for the six models that satisfy within one a all the observational 
constraints. The green squares show the observed values for Lx,LMXBs[2-l0ke\]/Mi, (top left panel) and Lx.HMXBs[Q.5-gkeV]/SFR (top right panel) as derived by 
Lehmeretal. (2010). The red squares show the observed constraints for Ax,LMA'fi.i[0.3-8keV]/A'* bv lBoroson et all (2011), and for £x.«MXBj[0.5-8keV]/SFK by 
Mineo et al. 1 2012a) after converted for demonstration purposes to the 2 - lOkeV energy band, in the left and right panel respectively. The bottom panel shows 
the X-ray luminosity density coming f rom the the whole XRB population as a function of redshift. The green squares show the observed values as reported by 
Tzanavaris & Georgantopoulos (2008), which we used as upper limits in our statistical analysis. The parameters of the different models are shown in Table[5] 

-,4 i ,_ 1 ____ r ___ r _ T __^____ r _, lldiart et all2007l for a discussion on this). 

In conclusion, we believe that further work will be needed, 
both observational in order to derive more reliable star- 
formation history information and theoretical incorporating 
dynamically formed LMXBs in our models, to unambigu- 
ously determine whether the sharp increase in Lx /M* going 
to ages younger than ~ 3 Gyr is real. This is clearly outside 
the scope of this study. However, the success of our models 
in predicting local Lx I M» and Lx I SFR re lations and the evo- 
lution of Lx /SFR (B asu-Zvch et al.ll2012l) while be consistent 
at the same time with X-ray stackin g results of the evo lution 
of Lx /M* from early-type galaxies dLehmer et al.H2007l) sug- 
gests that we are on the right track. 

5. RESULTS 

5.1. Effects of the different model parameters on the X-ray 
binary population 

As already mentioned, our models have two types of pa- 
rameters, the ones that correspond to distributions of initial 
properties of the stellar population, and which are guided by 
observations, and the "true" free parameters, that correspond 
to poorly understood physical processes which we are not able 
to model in detail. For both types of parameters, it is instruc- 
tive to study the effect they have in the XRB population and 
the constraints we can put on the values of these parameters 
by comparing the predictions of our models with observa- 




6 8 
Age (Gyr) 

Figure 6. Specific X-ray luminosity in the 0.3 — 8keV energy range, i.e. lu- 
minosity per unit stellar mass, of a coeval stellar population as a function of 
the population's age for our "reference" model (Model 245; see also sections 
2.4 and 4.2) and for solar metallicity. On the same figure over-plotted are 
the specific X-ray luminosities, in the same energy band, of nearby elliptical 
galaxies that have available in the literature estimates of their stellar popula- 
tion age. These data points are derived from the galaxy sample presented in 
Boros onetaUf20TTh . 

breaks down. This makes the ages highly uncertain, as a small 
fraction of young stars (< 10%) can make a galaxy with a 
much older population appear only a few Gyr old (see, e.g., 
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tions. For the first type of parameters, this comparison allows 
us to infer which part of the observationally allowed param- 
eter space is favored by our models. Constraining the sec- 
ond group of parameters is more essential, as we improve our 
understanding of complex physical processes that we cannot 
model in detail. 

Figure|7]shows the X-ray luminosity of the LMXB popula- 
tion, per unit stellar mass (left panel) and the X-ray luminos- 
ity of the HMXB population, per unit SFR (right panel), as 
a function of redshift, for eight selected models. The models 
shown in this figure were selected so that they demonstrate 
how each of the model parameters we varied in our grid af- 
fects the LMXB and HMXB population respectively. 

The CE efficiency a ce , i.e. the efficiency with which or- 
bital energy is converted into thermal energy that will expel 
the donor star's envelope during the CE phase, is a parame- 
ter that mainly affects the formation rate of LMXBs. Most 
LMXBs formed in the field, through the evolution of primor- 
dial binaries, have to go through a CE phase in order to form 
a tight enough binary that will later on become a LMXB. A 
higher a ce would result in more systems surviving the CE 
and forming tight post-CE binaries, and hence a higher for- 
mation rate of LMXBs. On the other hand, HMXBs do not 
necessarily have to go through a CE, but can instead become 
bright HMXBs throu gh alternative formation channels (e.g. 
IValsecchi etalfe OlO). Model 246 has the same model param- 
eters as our maximum likelihood model {Model 245) except 
for the higher CE efficiency. Comparing the two models (see 
Figure |7), we observe that indeed the LMXB population is 
heavily affected by a small change in a ce , while the HMXB 
population is rather insensitive. Our parameter study shows 
that only a low CE efficiency (a ce ~ 0.1) produces XRB pop- 
ulations consistent with observations. This finding may alter- 
natively be interpreted in terms of the binding energy param- 
eter A. Instead of the assumed value (A = 1) this parameter 
may be much lower (lambda ^0.1, if a ce is assumed 1) in- 
dicating that at the onset of the CE phase LMXB progenitors 
have highly bound envelopes. This in fact seems to be true 
for massive stars based on the actual estimates of e nvelope 
binding energies (see Fig. 1-4 of Domin ik et al.ll2012l) . 

Changing the initial binary mass ratio distribution between 
a flat distribution and 50% flat - 50% twin distribution affects 
the XRB population in a similar way (compare Model 53 and 
Model 245 in Figure[7J. The progenitor of a LMXB is a binary 
with high mass ratio. By forcing half of the systems to have 
an initial mass ratio close to 1, the formation rate of LMXBs 
is decreased by a factor of two. On the other hand, the HMXB 
population is largely unaffected. This suggests that in a binary 
population with a flat initial mass ratio distribution, approxi- 
mately 10% of HMXBs have as progenitors primordial bina- 
ries with a mass ratio between 0.9 and 1. The statistical com- 
parison of our models with observations favors a mixed initial 
mass ratio distribution of 50% flat - 50% twin distribution, 
which is consistent with observational surveys of young bina- 
ries (e.g. iPinsonneault &~S tanek 2006), although some recent 
studie s hint towards a single flat distribution dSana & Evansl 
1201 lh . A pure twin mass ratio distribution, although it is in- 
cluded in our parameter study (Models 100-199), results in 
models that are inconsistent with observations, as it will com- 
pletely prevent the formation of LMXBs. Figure |4]shows that 
models with a pure twin mass ratio distribution (Models 100- 
199) have systematically lower likelihood values than the rest. 

Stellar winds play also an important role in the formation 
and evolution of XRBs, and they do so primarily in two ways. 



A stronger stellar wind will increase the accretion rate on the 
compact object of a wind-fed HMXB, making it more lumi- 
nous. On the contrary, a weaker stellar wind will result in 
smaller overall mass loss of the primary star that will form 
the compact object, and hence in the formation of more nu- 
merous and more massive BHs. BH-XRBs tend to be more 
luminous than XRBs with NS accretors, as on the one hand 
they can form stable RLO XRB systems with massive com- 
panion stars, and on the other hand they can drive in gen- 
eral higher accretion rates (compared again to NS accretors), 
both in RLO and wind-fed systems, because of the higher 
mass of the BH accretor. Furthermore, a reduced stellar wind 
would result in less orbital expansion due to angular momen- 
tum loss from the wind, and in overall more systems that will 
encounter Roche-lobe overflow mass transfer. It turns out 
that the latter two effects are the dominant ones. Compar- 
ing Model 253 (rjwind = 2.0) and Model 261 (r/wind = 0.25) to 
our maximum likelihood model (Model 245, riwind = 1 -0), one 
can see the large effects: for stronger winds (model 253) the 
LMXB emission and HMXB emission is suppressed and for 
weaker winds (model 261) both the LMXB and HMXB emis- 
sion is boosted. Our statistical analysis showed that the "stan- 
da rd" wind prescript i on {riw ind = 1 0), compiled and calibrated 
by Be lczvnski et al. (2010), is favored, but models with stel- 
lar winds increased by a factor of two are also consistent with 
observations. If this result is taken at the face value (i.e., if po- 
tential degeneracies with other free parameters are neglected) 
it indicates that there is not much more unaccounted clumping 
in winds of massive stars that in the past led to overestimates 
of mass loss rates. 

The power law slope at the high-mass end of the IMF af- 
fects the population of XRBs in a way equivalent to stellar 
winds, in the sense that a flatter IMF will produce relatively 
more BHs compared to a steeper one. Hence, a flatter IMF 
will result in more luminous populations of both LMXBs and 
HMXBs (compare Model 241 to Model 245). In our standard 
treatment of supernova events, the most massive BHs receive 
no kick during their formation. Allowing for a small kick 
(10% of the standard kicks that NSs receive) in these systems 
does not affect significantly the LMXB population, but re- 
sults in an overall less luminous HMXB population (compare 
Model 245 to Model 197). HMXBs originate from more mas- 
sive binaries on average, where the formation of a BH through 
direct core-collapse is more likely. Thus, allowing for small 
kicks in these supernova events leads to the disruption of bi- 
naries that would otherwise form and HMXB. 

Finally, we found that allowing or not for all possible out- 
comes in CE phases with donor stars in the Hertzsprung gap, 
although it affec ts the shape of the XLF for a given galaxy 
dLuo et al.ll2012l) . has a negligible effect on the integrated X- 
ray luminosity of a LMXB population and slightly increases 
the luminosity of a HMXB population (compare Model 245 
to Model 269). This comes from the fact that LMXB pro- 
genitors are on average less massive than stars producing 
HMXBs. The probability of a CE phase to happen while 
the donor star is on the Hertzsprung gap is very low, since 
only very massive stars experience significant expansion dur- 
ing Hertzsprung gap. In contrast, for massive binaries the 
Hertzsprung gap CE happens frequently and if the survival 
is allowed some systems survive and form extra population of 
HMXBs. This additional sources lead to a factor of ~ 1.5 in- 
crease in the integrated HMXB X-ray luminosity throughout 
the cosmic time. This is very similar effect as noted for dou- 
ble compact object populations, for which massive BH-BH 
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Figure 7. Bolometric X-ray luminosity of the LMXB population, per unit stellar mass (left panel) and X-ray luminosity of the HMXB population, per unit SFR 
(right panel), as a function of redshift, for eight selected models that show the dependence of the models' predictions to the different modelparameters. Our 
maximum likelihood model (Model 245) is shown with the bold solid (black) line. The parameters of the different models are shown in Table|3] 



progenitors are greatly affected by Hertzsprung gap CE wh ile 
lighter NS-NS progenitors are not dBelczvnski et al.ll2007l) . 

5.2. Evolution of global scaling relations with redshift 

Figure [8] shows the bolometric X-ray luminosity per unit 
volume, unit of stellar mass, and unit of SFR, as a function 
of redshift, for our maximum likelihood model (Model 245). 
The relative contribution in these quantities of the LMXB and 
HMXB populations is also shown. Comparing in the top 
panel of Figure [8] the redshift evolution of the SFR density 
to the evolution of the X-ray luminosity density coming from 
LMXBs, one notices an evident time difference between the 
peaks of the two quantities. The peak of the SFR density oc- 
curs at z ~ 3.1, while the peak of the X-ray luminosity density 
coming from the LMXB populations happens at z ~ 2.1. This 
delay between the two peaks translates into a typical delay 
time of ~ 1.1 Gyr between the formation of a progenitor bi- 
nary star and time th at the systems reaches t he LMXB phase. 
In ear lier work by IWhite & Ghoshl (119981) : iGhosh & White! 
d2001l) . this timescale was a free parameter in their models 
which they we not able to constrain. We should note here, 
that the value of timescale is dependent on the exact defini- 
tion of what one calls a LMXB. As mentioned earlier, in this 
work we define as LMXB any RLO XRB with a donor star 
less massive than 3 M©. Shifting this mass boundary to lower 
masses, will lengthen the delay. 

Comparing now the SFR density peak to the peak of the 
X-ray luminosity density coming from HMXBs, one would 
expect that two peaks occur simultaneously. However, the 
peak X-ray luminosity density coming from HMXBs occurs 
at z ~ 3 .9, ~ 0.5 Gyr before the peak of the star formation. As 
we can see from the bottom panel of Figure [8] the X-ray lu- 
minosity coming from HMXBs, per unit SFR, is not constant 
with time, instead it is increasing with redshift. This is a re- 
sult of the metallicity evolution. The mean metallicity of the 
newly formed stars is decreasing with redshift (see Figure [T] 
lower panel). A lower metallicity translates into weaker stellar 
winds, which in turn results in more numerous and more mas- 
sive BHs. Since BH-HMXBs tend to be more luminous than 
NS-HMXBs, the X-ray luminosity coming from HMXBs per 
unit SFR is higher for lower metallicity stellar populations, 
and hence at higher redshifts. Coming back to the time differ- 
ence between the peak of the SFR and the peak of the X-ray 
luminosity density from HMXBs, it is now evident that it is 
because of the dependence of the HMXB population to the 



metallicity and the evolution of the metallicity with redshift. 

One additional point that is demonstrated in Figure [8] is 
that LMXBs dominate the X-ray emission from XRBs in our 
Universe today by a factor of ~ 2.5. The transition between 
an HMXB and LMXB dominated XRB population happens 
at a z ~ 2.5. This finding implies that the X-ray luminos- 
ity of an XRB population in a "typical" galaxy, with a star- 
formation history similar to the global star-formation his- 
tory of the Universe, is dominated by the X-ray emission of 
LMXBs. This is not an unexpected result, given that in the 
Milky Way, which is a typical spiral galaxy, ther e are only a 
few HMXBs with luminos ity above 10 37 ergs _1 ( Grim m et al.1 
l2002t IVoss & Ajellol 120101) . However, our finding should 
serve as an extra caution to a usual misconception that the 
XRB population of galaxies that currently sh ow some star 
forma tion have a negligible LMXB population. Lehm eret al.l 
(2010) showed observational evidence that the X-ray lumi- 
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nosity from galaxies with SFR as high as 
non-negligible contribution from LMXBs. 

Our modeling predicts that the X-ray luminosity coming 
from HMXBs, per unit SFR is increasing with redshift, and 
that this is a result of the metallicity evolution of the Uni- 
verse. The evolution of the X-ray luminosity from HMXBs 
per unit SFR predicted by our models (by a factor of ~ 10 
out to z ~ 15) is consistent with constrain t set by th e ob- 
served cosmic X-ray background (Diik stra et al.ll2012l) . Al- 
though the X-ray luminosity coming from HMXBs, per unit 
SFR, changes by a factor of 5 between z = and z ~ 10, the 
observed signature of this effect would be masked by the con- 
tamination of LMXBs. The bottom panel of Figure [8] shows 
that the contamination of LMXBs operates in such a way that 
the X-ray luminosity from the whole XRB population, per 
unit SFR, remains approximately constant with t ime. These 
results are in agreement with t he findings by Lehm eret al.l 
(2008) and lMineo et al.l d2012bl) for the late- type galaxy pop- 
ulation. Hence, our modeling suggests that the effect of the 
metallicity evolution to the HMXB population can only be 
observed with a carefully selected galaxy sample, where all 
galaxies have high enough SFR to stellar mass ratios, ensur- 
ing that their XRB populations are dominated by HMXBs. 

The middle panel of Figure 7 shows that the X-ray lumi- 
nosity of the whole XRB population, per unit stellar mass, 
evolves with redshift following approximately a single power 
law from z = 0, back to z ~ 15. This relation probes two dis- 
tinct phases of the cosmic evolution of the XRB population. 
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At z < 2.5, where the emission from LMXBs dominates the 
total X-ray luminosity, the evolution of the X-ray luminos- 
ity per unit stellar mass is related to the average age of the 
LMXB population. The fact that relatively younger stellar 
populations are associated with more luminous populations of 
LMXBs has been already observed in X-ray stacking studies 
of distant galaxies dLehmer et al.ll2007l) . iLehmer et all (120071) 
investigated the average X-ray properties of early-type galax- 
ies within the Extended Chandra Deep Field-South, and found 
a slight increase of Lx/Lb,q, where Lb i0 is the optical lumi- 
nosity in the B-band, with redshift, over the redshift range 
Z = 0-0.7, in excellent agreement with our model predictions. 
At z > 2.5, it is the emission from HMXBs that dominates 
the X-ray luminosity of the XRB population. The total X- 
ray luminosity coming from HMXBs is to first order, neglect- 
ing the metallicity evolution effects, proportional to the SFR. 
Hence, the evolution of X-ray luminosity per unit stellar mass 
at z > 2.5, can be a probe of the evolution of the specific SFR 
(SFR per unit stellar mass) with redshift. 

In the interpretation of our simulation results, and the com- 
parison with the observational data, one should keep in mind 
that all the quantities calculated in this paper correspond to 
the average of all galaxies in the Universe. A more de- 
tailed comparison of our models with galaxy surveys of early 
and late type galaxies separately, where all the sample se- 
lection effects are taken int o account, is discussed in other 
works (iTremmel et al.ll2012t iBasu-Zvch elaTI 120121 ; Horn- 
schemeier et al. 2012, in prep.). 

6. SUMMARY AND CONCLUSIONS 

We present here the first PS study of XRBs in a cosmolog- 
ical context. We created the largest grid of 288 PS models 
over nine metallicities, where we varied all the major model 
parameters that are known to affect the evolution of XRBs. 
We convolved these models with information about the star- 
formation history and metallicity evolutio n of the Universe 
from the Millennium II simulation and the iGuo et all (1201 ll) 
semi-analytic galaxy catalog. The combination of our PS 
models with a cosmological simulation allowed us to derive 
global scaling relations of the XRB population with proper- 
ties such as SFR and stellar mass, and the evolution of these 
relations with redshift. These predictions were compared sta- 
tistically with observations from the local Universe in order to 
constrain our models and make robust predictions about the 
high-redshift Universe. The main conclusions of this work 
can be summarized as follows: 

1 . The statistical comparison of our model predictions re- 
vealed several models that are in excellent agreement 
with observations of the local Universe. Furthermore, 
we found that all models that are consistent with obser- 
vational data in the local Universe exhibit very similar 
behavior at higher redshifts, indicating that our predic- 
tions for the redshift evolution of XRB populations are 
robust. 

2. Although the main goal of this paper is not to constrain 
the parameters of our PS models but instead to study the 
redshift evolution of XRB populations, our parameter 
study allows for some general conclusions on the possi- 
ble values of the PS parameters we varied. The statisti- 
cal analysis showed that a CE efficiency A x a C£ < 0.1, 
where A is a measure of the central concentration of 
the donor star, and possibly a mixed initial mass ratio 
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Figure 8. Bolometric X-ray luminosity per unit volume (top panel), unit of 
stellar mass (middle panel), and unit of SFR (bottom panel), as a function of 
redshift, for our maximum likelihood model (Model 245). The solid (black) 
lines correspond to the whole XRB population, while the dotted (blue) lines 
show the contribution of the HMXBs, and the dashed (red) lines the contri- 
bution of the LMXBs. For comparison, the SFR density is shown on the top 
panel with a grey dash-dotted line. The vertical solid (black), dotted (blue), 
dashed (red), and dash-dotted (grey) lines denote the redshifts at which the 
peaks of the X-ray luminosity density from the whole XRB population, the 
HMXBs, and the LMXBs, and the SFR density occur. 

distribution (i.e. 50% of binaries are drawn from a flat 
initial mass ratio distribution and 50% from a '"twin" 
distributuon) is necessary in order to produce a LMXB 
population consistent with observations. These two pa- 
rameters show some degeneracy, and a slightly lower 
CE efficiency would possibly favor a flat initial mass 
ratio distribution. On the other hand, model parame- 
ters that regulate the formation rate of BH-XRBs, such 
as the IMF and the stellar wind strength, appear to be 
important for both the LMXB and HMXB populations. 
This is not surprising as XRBs with BH accretors can 
be very luminous, and even a relatively small number 
of them can dominate the integrated X-ray luminosity 
of a population. Again, these parameters show a de- 
generacy as they affect the population in a similar way, 
and only their combined effect can be constrained from 
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observations. 

3. The X-ray luminosity of XRBs per unit stellar mass or 
unit star-formation rate, is significantly higher for low- 
metallicity stellar populations, due to the fact that the 
young massive O/B stars, which are the precursors to 
compact object accretors in X-ray binaries, will lose 
less stellar mass from line-driven winds over their life- 
times for lower metallicities. This results in more nu- 
merous and more massive BHs, and subsequently more 
luminous X-ray binary populations. 

4. The X-ray emission from XRBs in the local Universe 
is dominated by LMXBs, and it is only at z > 2.5 that 
HMXBs start to dominate the X-ray emission. The red- 
shift that this transition happens is lower than that of the 
peak of star formation (z ~ 3), which shows that there 
is a delay between the formation of primordial binaries 
and the time at which these binaries become LMXBs. 
The formation timescale of LMXBs is estimated to be 
~ 1 . 1 Gyr. This finding also implies that the X-ray lu- 
minosity of an XRB population in a "typical" galaxy, 
with a star-formation history similar to the global star- 
formation history of the Universe, is dominated by the 
X-ray emission of LMXBs. 

5. The formation rate of HMXBs follows closely the SFR. 
However, the X-ray luminosity from HMXBs per unit 
SFR (Lx, hmxbs/ SFR) depends also on metallicity. Tak- 
ing into account the metallicity evolution of the Uni- 
verse, we find that the peak of the X-ray luminosity 
from HMXBs per unit volume, occurs ~ 0.5 Gyr before 
the peak of the SFR density. Although the X-ray lu- 
minosity coming from HMXBs, per unit SFR, changes 
by half an order of magnitude over the evolution of 
the Universe, the observed signature of this effect is 
masked by the contamination of LMXBs, making the 
X-ray luminosity from the whole XRB population per 
unit SFR approximately constant with redshift. 

6. Finally, we found that at z < 2.5, where the emission 
from LMXBs dominates the total X-ray luminosity, the 
evolution of the X-ray luminosity per unit stellar mass 
is related to the average age of the LMXB population, 
while at z > 2.5, the evolution of X-ray luminosity per 
unit stellar mass can be a probe of the evolution of the 
specific SFR with redshift. 

This first PS study of the evolution of XRB populations 
across comic time highlights the importance of the study of X- 
ray luminous normal galaxies at high redshift. Furthermore, 
it lays the ground for future work that will study role played 
by XRBs in the formatio n and evolution of galaxies through 
feedback processes (e. g. Mirabel et al.l 120111: iHaimanl 1201 It 
Uustham & Schawinskill2012l) . 
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APPENDIX 

TO APPEAR AS SUPPLEMENTAL ONLINE-ONLY MATERIAL 



T. Fragos et al. 



Table 4 

A complete list of all PS models used in this work. 



Model 




IMF exponent 




CE-HG d 


q distribution 6 


K-DCBH f 


Ranke 


log(Likelihood ratio) h 


1 


0.1 


-2.35 


1.0 


No 


flat 


0.0 


81 


-6.8439794 


2 


0.2 


-2.35 


1.0 


No 


flat 


0.0 


131 


-13.449294 


3 


0.3 


-2.35 


1.0 


No 


flat 


0.0 


149 


-17.435949 


4 


0.5 


-2.35 


1.0 


No 


flat 


0.0 


167 


-22.784369 


5 


0.1 


-2.7 


1.0 


No 


flat 


0.0 


23 


-1.7703540 


6 


0.2 


-2.7 


1.0 


No 


flat 


0.0 


65 


-5.9952885 


7 


0.3 


-2.7 


1.0 


No 


flat 


0.0 


99 


-8.6721399 


8 


0.5 


-2.7 


1.0 


No 


flat 


0.0 


121 


-11.957086 


9 


0.1 


-2.35 


2.0 


No 


flat 


0.0 


29 


-2.3030907 


10 


0.2 


-2.35 


2.0 


No 


flat 


0.0 


93 


-8.1440863 


11 


0.3 


-2.35 


2.0 


No 


flat 


0.0 


113 


-11.281685 


12 


0.5 


-2.35 


2.0 


No 


flat 


0.0 


115 


-11.393684 


13 


0.1 


-2.7 


2.0 


No 


flat 


0.0 


9 


-0.50583441 


14 


0.2 


-2.7 


2.0 


No 


flat 


0.0 


51 


-4.3697639 


15 


0.3 


-2.7 


2.0 


No 


flat 


0.0 


70 


-6.2810696 


16 


0.5 


-2.7 


2.0 


No 


flat 


0.0 


64 


-5.8452086 


17 


0.1 


-2.35 


0.25 


No 


flat 


0.0 


160 


-20.757946 


18 


0.2 


-2.35 


0.25 


No 


flat 


0.0 


175 


-27.096760 


19 


0.3 


-2.35 


0.25 


No 


flat 


0.0 


183 


-32.066400 


20 


0.5 


-2.35 


0.25 


No 


flat 


0.0 


190 


-40.521407 


21 


0.1 


-2.7 


0.25 


No 


flat 


0.0 


101 


-8.9727405 


22 


0.2 


-2.7 


0.25 


No 


flat 


0.0 


137 


-14.797751 


23 


0.3 


-2.7 


0.25 


No 


flat 


0.0 


151 


-18.195038 


24 


0.5 


-2.7 


0.25 


No 


flat 


0.0 


168 


-23.306565 


25 


0.1 


-2.35 


1.0 


Yes 


flat 


0.0 


94 


-8.1863846 


26 


0.2 


-2.35 


1.0 


Yes 


flat 


0.0 


138 


-14.962405 


27 


0.3 


-2.35 


1.0 


Yes 


flat 


0.0 


156 


-19.167849 


28 


0.5 


-2.35 


1.0 


Yes 


flat 


0.0 


174 


-26.883317 


29 


0.1 


-2.7 


1.0 


Yes 


flat 


0.0 


24 


-1.8904716 


30 


0.2 


-2.7 


1.0 


Yes 


flat 


0.0 


82 


-6.8828256 


31 


0.3 


-2.7 


1.0 


Yes 


flat 


0.0 


108 


-10.076327 


32 


0.5 


-2.7 


1.0 


Yes 


flat 


0.0 


134 


-14.094516 


33 


0.1 


-2.35 


2.0 


Yes 


flat 


0.0 


30 


-2.3061373 


34 


0.2 


-2.35 


2.0 


Yes 


flat 


0.0 


100 


-8.8750378 


35 


0.3 


-2.35 


2.0 


Yes 


flat 


0.0 


122 


-12.034528 


36 


0.5 


-2.35 


2.0 


Yes 


flat 


0.0 


128 


-13.072784 


37 


0.1 


-2.7 


2.0 


Yes 


flat 


0.0 


7 


-0.45466309 


38 


0.2 


-2.7 


2.0 


Yes 


flat 


0.0 


49 


-4.1820059 


39 


0.3 


-2.7 


2.0 


Yes 


flat 


0.0 


76 


-6.4496818 


40 


0.5 


-2.7 


2.0 


Yes 


flat 


0.0 


74 


-6.4069133 


41 


0.1 


-2.35 


0.25 


Yes 


flat 


0.0 


159 


-20.659322 


42 


0.2 


-2.35 


0.25 


Yes 


flat 


0.0 


182 


-30.351627 


43 


0.3 


-2.35 


0.25 


Yes 


flat 


0.0 


187 


-36.893741 


44 


0.5 


-2.35 


0.25 


Yes 


flat 


0.0 


192 


-48.616001 


45 


0.1 


-2.7 


0.25 


Yes 


flat 


0.0 


106 


-9.8900460 


46 


0.2 


-2.7 


0.25 


Yes 


flat 


0.0 


143 


-16.465170 


47 


0.3 


-2.7 


0.25 


Yes 


flat 


0.0 


164 


-21.251989 


48 


0.5 


-2.7 


0.25 


Yes 


flat 


0.0 


179 


-28.630340 


49 


0.1 


-2.35 


1.0 


No 


flat 


0.1 


58 


-5.4944514 


50 


0.2 


-2.35 


1.0 


No 


flat 


0.1 


119 


-11.860356 


51 


0.3 


-2.35 


1.0 


No 


flat 


0.1 


140 


-15.604207 


52 


0.5 


-2.35 


1.0 


No 


flat 


0.1 


165 


-21.300323 


53 


0.1 


-2.7 


1.0 


No 


flat 


0.1 


20 


-1.3945094 


54 


0.2 


-2.7 


1.0 


No 


flat 


0.1 


59 


-5.5543889 


55 


0.3 


-2.7 


1.0 


No 


flat 


0.1 


96 


-8.3582206 


56 


0.5 


-2.7 


1.0 


No 


flat 


0.1 


109 


-10.902530 


57 


0.1 


-2.35 


2.0 


No 


flat 


0.1 


28 


-2.2146716 


58 


0.2 


-2.35 


2.0 


No 


flat 


0.1 


92 


-8.0877140 


59 


0.3 


-2.35 


2.0 


No 


flat 


0.1 


112 


-11.149264 


60 


0.5 


-2.35 


2.0 


No 


flat 


0.1 


111 


-11.047687 


61 


0.1 


-2.7 


2.0 


No 


flat 


0.1 


16 


-1.1282696 


62 


0.2 


-2.7 


2.0 


No 


flat 


0.1 


52 


-4.4089974 


63 


0.3 


-2.7 


2.0 


No 


flat 


0.1 


77 


-6.5627561 


64 


0.5 


-2.7 


2.0 


No 


flat 


0.1 


72 


-6.3313205 


OJ 


U. 1 


-Z.JJ 


VI. z. J 


1NO 


flat 

nat 


U. 1 


1/17 


- lo.y /Uo jZ 


66 


0.2 


-2.35 


0.25 


No 


flat 


0.1 


170 


-24.339411 


67 


0.3 


-2.35 


0.25 


No 


flat 


0.1 


181 


-28.872144 


68 


0.5 


-2.35 


0.25 


No 


flat 


0.1 


188 


-37.108159 


69 


0.1 


-2.7 


0.25 


No 


flat 


0.1 


79 


-6.6803483 


70 


0.2 


-2.7 


0.25 


No 


flat 


0.1 


123 


-12.306121 


71 


0.3 


-2.7 


0.25 


No 


flat 


0.1 


141 


-15.707806 


72 


0.5 


-2.7 


0.25 


No 


flat 


0.1 


163 


-21.235448 


73 


0.1 


-2.35 


1.0 


Yes 


flat 


0.1 


71 


-6.3013718 


74 


0.2 


-2.35 


1.0 


Yes 


flat 


0.1 


129 


-13.103548 


75 


0.3 


-2.35 


1.0 


Yes 


flat 


0.1 


153 


-18.401280 
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Table 4 — Continued 



Model 


a 

&CE 


IMr exponent 


c 

IJwind 


V [-[Hi 


r i ■ e 
q distribution 


f 

K-DCBH 


Pernio 
RdilK" 


. ... . h 
log^LiKelinood ratio) 


76 


0.5 


-2.35 


1.0 


Yes 


flat 


0.1 


172 


-25.574837 


77 


0.1 


-2.7 


1.0 


Yes 


flat 


0.1 


21 


-1.5355082 


78 


0.2 


-2.7 


1.0 


Yes 


flat 


0.1 


73 


-6.4029360 


79 


0.3 


-2.7 


1.0 


Yes 


flat 


0.1 


103 


-9.3432712 


80 


0.5 


-2.7 


1.0 


Yes 


flat 


0.1 


130 


-13.387239 


81 


0.1 


-2.35 


2.0 


Yes 


flat 


0.1 


26 


-2.1743665 


82 


0.2 


-2.35 


2.0 


Yes 


flat 


0.1 


98 


-8.5541878 


83 


0.3 


-2.35 


2.0 


Yes 


flat 


0.1 


118 


-11.762885 


84 


0.5 


-2.35 


2.0 


Yes 


flat 


0.1 


127 


-13.021688 


85 


0.1 


-2.7 


2.0 


Yes 


flat 


0.1 


12 


-0.78509022 


86 


0.2 


-2.7 


2.0 


Yes 


flat 


0.1 


50 


-4.2847568 


87 


0.3 


-2.7 


2.0 


Yes 


flat 


0.1 


68 


-6.2687203 


88 


0.5 


-2.7 


2.0 


Yes 


flat 


0.1 


75 


-6.4135880 


89 


0.1 


-2.35 


0.25 


Yes 


flat 


0.1 


152 


-18.205247 


90 


0.2 


-2.35 


0.25 


Yes 


flat 


0.1 


176 


-27.999553 


91 


0.3 


-2.35 


0.25 


Yes 


flat 


0.1 


186 


-34.997344 


92 


0.5 


-2.35 


0.25 


Yes 


flat 


0.1 


191 


-46.958394 


93 


0.1 


-2.7 


0.25 


Yes 


flat 


0.1 


86 


-7.1404740 


94 


0.2 


-2.7 


0.25 


Yes 


flat 


0.1 


136 


-14.372125 


95 


0.3 


-2.7 


0.25 


Yes 


flat 


0.1 


158 


-20.101154 


96 


0.5 


-2.7 


0.25 


Yes 


flat 


0.1 


178 


-28.160352 


97 


0.1 


-2.35 


1.0 


No 


twin 


0.0 


260 


-156.87865 


98 


0.2 


-2.35 


1.0 


No 


twin 


0.0 


226 


-118.90171 


99 


0.3 


-2.35 


1.0 


No 


twin 


0.0 


227 


-120.00800 


100 


0.5 


-2.35 


1.0 


No 


twin 


0.0 


273 


-173.98895 


101 


0.1 


-2.7 


1.0 


No 


twin 


0.0 


274 


-177.22909 


102 


0.2 


-2.7 


1.0 


No 


twin 


0.0 


234 


-126.11918 


103 


0.3 


-2.7 


1.0 


No 


twin 


0.0 


239 


-131.17535 


104 


0.5 


-2.7 


1.0 


No 


twin 


0.0 


280 


-189.56799 


105 


0.1 


-2.35 


2.0 


No 


twin 


0.0 


276 


-181.10990 


106 


0.2 


-2.35 


2.0 


No 


twin 


0.0 


246 


-139.37974 


107 


0.3 


-2.35 


2.0 


No 


twin 


0.0 


232 


-124.52035 


108 


0.5 


-2.35 


2.0 


No 


twin 


0.0 


279 


-183.77797 


109 


0.1 


-2.7 


2.0 


No 


twin 


0.0 


283 


-212.54419 


110 


0.2 


-2.7 


2.0 


No 


twin 


0.0 


249 


-141.81642 


111 


0.3 


-2.7 


2.0 


No 


twin 


0.0 


252 


-144.96094 


112 


0.5 


-2.7 


2.0 


No 


twin 


0.0 


275 


-180.61817 


113 


0.1 


-2.35 


0.25 


No 


twin 


0.0 


281 


-190.13124 


114 


0.2 


-2.35 


0.25 


No 


twin 


0.0 


209 


-105.58207 


115 


0.3 


-2.35 


0.25 


No 


twin 


0.0 


220 


-113.36827 


116 


0.5 


-2.35 


0.25 


No 


twin 


0.0 


243 


-136.34902 


117 


0.1 


-2.7 


0.25 


No 


twin 


0.0 


271 


-173.34531 


118 


0.2 


-2.7 


0.25 


No 


twin 


0.0 


212 


-107.65630 


119 


0.3 


-2.7 


0.25 


No 


twin 


0.0 


223 


-115.87870 


120 


0.5 


-2.7 


0.25 


No 


twin 


0.0 


259 


-156.69007 


121 


0.1 


-2.35 


1.0 


Yes 


twin 


0.0 


265 


-166.58676 


122 


0.2 


-2.35 


1.0 


Yes 


twin 


0.0 


236 


-128.34509 


123 


0.3 


-2.35 


1.0 


Yes 


twin 


0.0 


202 


-102.01665 


124 


0.5 


-2.35 


1.0 


Yes 


twin 


0.0 


199 


-99.689704 


125 


0.1 


-2.7 


1.0 


Yes 


twin 


0.0 


270 


-171.79474 


126 


0.2 


-2.7 


1.0 


Yes 


twin 


0.0 


245 


-138.27339 


127 


0.3 


-2.7 


1.0 


Yes 


twin 


0.0 


231 


-123.30591 


128 


0.5 


-2.7 


1.0 


Yes 


twin 


0.0 


219 


-111.26572 


129 


0.1 


-2.35 


2.0 


Yes 


twin 


0.0 


269 


-170.91556 


130 


0.2 


-2.35 


2.0 


Yes 


twin 


0.0 


244 


-138.08749 


131 


0.3 


-2.35 


2.0 


Yes 


twin 


0.0 


233 


-125.38742 


132 


0.5 


-2.35 


2.0 


Yes 


twin 


0.0 


204 


-103.05276 


133 


0.1 


-2.7 


2.0 


Yes 


twin 


0.0 


287 


-238.92812 


134 


0.2 


-2.7 


2.0 


Yes 


twin 


0.0 


253 


-145.13205 


135 


0.3 


-2.7 


2.0 


Yes 


twin 


0.0 


250 


-142.40690 


136 


0.5 


-2.7 


2.0 


Yes 


twin 


0.0 


225 


-117.78790 


137 


0.1 


-2.35 


0.25 


Yes 


twin 


0.0 


257 


-150.66356 


138 


0.2 


-2.35 


0.25 


Yes 


twin 


0.0 


206 


-103.42077 


139 


0.3 


-2.35 


0.25 


Yes 


twin 


0.0 


197 


-99.321415 


140 


0.5 


-2.35 


0.25 


Yes 


twin 


0.0 


198 


-99.520795 


141 


0.1 


-2.7 


0.25 


Yes 


twin 


0.0 


268 


-170.03313 


142 


0.2 


-2.7 


0.25 


Yes 


twin 


0.0 


221 


-115.28205 


143 


0.3 


-2.7 


0.25 


Yes 


twin 


0.0 


196 


-99.065831 


144 


0.5 


-2.7 


0.25 


Yes 


twin 


0.0 


203 


-103.00718 


145 


0.1 


-2.35 


1.0 


No 


twin 


0.1 


282 


-212.05569 


146 


0.2 


-2.35 


1.0 


No 


twin 


0.1 


211 


-107.37148 


147 


0.3 


-2.35 


1.0 


No 


twin 


0.1 


215 


-109.36328 


148 


0.5 


-2.35 


1.0 


No 


twin 


0.1 


247 


-140.54672 


149 


0.1 


-2.7 


1.0 


No 


twin 


0.1 


261 


-158.65389 


150 


0.2 


-2.7 


1.0 


No 


twin 


0.1 


240 


-131.79837 


151 


0.3 


-2.7 


1.0 


No 


twin 


0.1 


228 


-120.63145 


152 


0.5 


-2.7 


1.0 


No 


twin 


0.1 


263 


-163.96495 



T. Fragos et al. 



Table 4 — Continued 



Model 


a 

&CE 


IMr exponent 


c 

Vwind 


pp i_rr;d 

V 1 -[ R 1 


r i ■ e 
q distribution 


f 

KDCBH 


Pernio 
RdilK" 


. ... . h 
log(L.ikelinoocl ratio) 


153 


0.1 


-2.35 


2.0 


No 


twin 


0.1 


285 


-226.73511 


154 


0.2 


-2.35 


2.0 


No 


twin 


0.1 


255 


-147.27372 


155 


0.3 


-2.35 


2.0 


No 


twin 


0.1 


241 


-131.82939 


156 


0.5 


-2.35 


2.0 


No 


twin 


0.1 


278 


-183.64930 


157 


0.1 


-2.7 


2.0 


No 


twin 


0.1 


286 


-229.31936 


158 


0.2 


-2.7 


2.0 


No 


twin 


0.1 


272 


-173.45129 


159 


0.3 


-2.7 


2.0 


No 


twin 


0.1 


256 


-149.25108 


160 


0.5 


-2.7 


2.0 


No 


twin 


0.1 


284 


-213.83915 


161 


0.1 


-2.35 


0.25 


No 


twin 


0.1 


267 


-169.73810 


162 


0.2 


-2.35 


0.25 


No 


twin 


0.1 


200 


-101.25795 


163 


0.3 


-2.35 


0.25 


No 


twin 


0.1 


216 


-109.61989 


164 


0.5 


-2.35 


0.25 


No 


twin 


0.1 


237 


-128.62607 


165 


0.1 


-2.7 


0.25 


No 


twin 


0.1 


254 


-146.26185 


166 


0.2 


-2.7 


0.25 


No 


twin 


0.1 


210 


-105.95448 


167 


0.3 


-2.7 


0.25 


No 


twin 


0.1 


235 


-127.40935 


168 


0.5 


-2.7 


0.25 


No 


twin 


0.1 


251 


-142.65248 


169 


0.1 


-2.35 


1.0 


Yes 


twin 


0.1 


277 


-182.07044 


170 


0.2 


-2.35 


1.0 


Yes 


twin 


0.1 


224 


-116.16233 


171 


0.3 


-2.35 


1.0 


Yes 


twin 


0.1 


217 


-110.39158 


172 


0.5 


-2.35 


1.0 


Yes 


twin 


0.1 


201 


-101.26795 


173 


0.1 


-2.7 


1.0 


Yes 


twin 


0.1 


213 


-108.61106 


174 


0.2 


-2.7 


1.0 


Yes 


twin 


0.1 


230 


-122.80642 


175 


0.3 


-2.7 


1.0 


Yes 


twin 


0.1 


222 


-115.68776 


176 


0.5 


-2.7 


1.0 


Yes 


twin 


0.1 


218 


-111.09642 


177 


0.1 


-2.35 


2.0 


Yes 


twin 


0.1 


262 


-158.68953 


178 


0.2 


-2.35 


2.0 


Yes 


twin 


0.1 


248 


-141.13179 


179 


0.3 


-2.35 


2.0 


Yes 


twin 


0.1 


238 


-130.61604 


180 


0.5 


-2.35 


2.0 


Yes 


twin 


0.1 


207 


-103.52411 


181 


0.1 


-2.7 


2.0 


Yes 


twin 


0.1 


288 


-254.52814 


182 


0.2 


-2.7 


2.0 


Yes 


twin 


0.1 


264 


-166.30555 


183 


0.3 


-2.7 


2.0 


Yes 


twin 


0.1 


242 


-132.91527 


184 


0.5 


-2.7 


2.0 


Yes 


twin 


0.1 


229 


-121.68847 


185 


0.1 


-2.35 


0.25 


Yes 


twin 


0.1 


266 


-169.59003 


186 


0.2 


-2.35 


0.25 


Yes 


twin 


0.1 


208 


-103.86318 


187 


0.3 


-2.35 


0.25 


Yes 


twin 


0.1 


194 


-94.819491 


188 


0.5 


-2.35 


0.25 


Yes 


twin 


0.1 


193 


-92.661723 


189 


0.1 


-2.7 


0.25 


Yes 


twin 


0.1 


258 


-156.51391 


190 


0.2 


-2.7 


0.25 


Yes 


twin 


0.1 


214 


-108.86360 


191 


0.3 


-2.7 


0.25 


Yes 


twin 


0.1 


205 


-103.13924 


192 


0.5 


-2.7 


0.25 


Yes 


twin 


0.1 


195 


-98.010507 


193 


0.1 


-2.35 


1.0 


No 


50-50 


0.0 


56 


-4.6196089 


194 


0.2 


-2.35 


1.0 


No 


50-50 


0.0 


102 


-8.9811857 


195 


0.3 


-2.35 


1.0 


No 


50-50 


0.0 


120 


-11.949863 


196 


0.5 


-2.35 


1.0 


No 


50-50 


0.0 


142 


-16.147741 


197 


0.1 


-2.7 


1.0 


No 


50-50 


0.0 


13 


-0.80759106 


198 


0.2 


-2.7 


1.0 


No 


50-50 


0.0 


39 


-2.8435160 


199 


0.3 


-2.7 


1.0 


No 


50-50 


0.0 


55 


-4.5791773 


200 


0.5 


-2.7 


1.0 


No 


50-50 


0.0 


83 


-6.9038795 


201 


0.1 


-2.35 


2.0 


No 


50-50 


0.0 


10 


-0.59598324 


202 


0.2 


-2.35 


2.0 


No 


50-50 


0.0 


48 


-3.9985781 


203 


0.3 


-2.35 


2.0 


No 


50-50 


0.0 


66 


-6.1997086 


204 


0.5 


-2.35 


2.0 


No 


50-50 


0.0 


69 


-6.2738874 


205 


0.1 


-2.7 


2.0 


No 


50-50 


0.0 


4 


-0.13235138 


206 


0.2 


-2.7 


2.0 


No 


50-50 


0.0 


18 


-1.3037351 


207 


0.3 


-2.7 


2.0 


No 


50-50 


0.0 


32 


-2.4529087 


208 


0.5 


-2.7 


2.0 


No 


50-50 


0.0 


27 


-2.1765207 


209 


0.1 


-2.35 


0.25 


No 


50-50 


0.0 


145 


-16.868460 


210 


0.2 


-2.35 


0.25 


No 


50-50 


0.0 


166 


-21.835832 


211 


0.3 


-2.35 


0.25 


No 


50-50 


0.0 


173 


-25.882512 


212 


0.5 


-2.35 


0.25 


No 


50-50 


0.0 


184 


-32.902336 


213 


0.1 


-2.7 


0.25 


No 


50-50 


0.0 


84 


-6.9688649 


214 


0.2 


-2.7 


0.25 


No 


50-50 


0.0 


110 


-10.931015 


215 


0.3 


-2.7 


0.25 


No 


50-50 


0.0 


132 


-13.487211 


216 


0.5 


-2.7 


0.25 


No 


50-50 


0.0 


150 


-17.481734 


217 


0.1 


-2.35 


1.0 


Yes 


50-50 


0.0 


60 


-5.5770783 


218 


0.2 


-2.35 


1.0 


Yes 


50-50 


0.0 


107 


-9.9479274 


219 


0.3 


-2.35 


1.0 


Yes 


50-50 


0.0 


126 


-12.981021 


220 


0.5 


-2.35 


1.0 


Yes 


50-50 


0.0 


154 


-18.960596 


221 


0.1 


-2.7 


1.0 


Yes 


50-50 


0.0 


15 


-1.0050246 


222 


0.2 


-2.7 


1.0 


Yes 


50-50 


0.0 


42 


-3.4137746 


223 


0.3 


-2.7 


1.0 


Yes 


50-50 


0.0 


57 


-5.4813514 


224 


0.5 


-2.7 


1.0 


Yes 


50-50 


0.0 


95 


-8.3259810 


225 


0.1 


-2.35 


2.0 


Yes 


50-50 


0.0 


11 


-0.75550478 


226 


0.2 


-2.35 


2.0 


Yes 


50-50 


0.0 


54 


-4.5531285 


227 


0.3 


-2.35 


2.0 


Yes 


50-50 


0.0 


80 


-6.7609447 


228 


0.5 


-2.35 


2.0 


Yes 


50-50 


0.0 


89 


-7.4758766 


229 


0.1 


-2.7 


2.0 


Yes 


50-50 


0.0 


2 


-0.035048515 
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IVlUUcl 


acE 


Tl\/TTn ovnnn on 
liVLF CApUIlCIlL 




CE-HG^ 


c| distribution 


f 

KDCBH 


Rank g 


i n i YU j f- \h 

10g(ljlKClinOOQ IcltlO) 


230 


0.2 


-2.7 


2.0 


Yes 


50-50 


0.0 


17 


-1.2284989 


231 


0.3 


-2.7 


2.0 


Yes 


50-50 


0.0 


37 


-2.6011593 


232 


0.5 


-2.7 


2.0 


Yes 


50-50 


0.0 


35 


-2.5733454 


233 


0.1 


-2.35 


0.25 


Yes 


50-50 


0.0 


144 


-16.735734 


234 


0.2 


-2.35 


0.25 


Yes 


50-50 


0.0 


169 


-23.741728 


235 


0.3 


-2.35 


0.25 


Yes 


50-50 


0.0 


180 


-28.856645 


236 


0.5 


-2.35 


0.25 


Yes 


50-50 


0.0 


189 


-38.542077 


237 


0.1 


-2.7 


0.25 


Yes 


50-50 


0.0 


91 


-7.5587255 


238 


0.2 


-2.7 


0.25 


Yes 


50-50 


0.0 


117 


-11.759989 


239 


0.3 


-2.7 


0.25 


Yes 


50-50 


0.0 


139 


-15.228415 


240 


0.5 


-2.7 


0.25 


Yes 


50-50 


0.0 


161 


-20.960307 


241 


0.1 


-2.35 


1.0 


No 


50-50 


0.1 


31 


-2.3582808 


242 


0.2 


-2.35 


1.0 


No 


50-50 


0.1 


78 


-6.6277666 


243 


0.3 


-2.35 


1.0 


No 


50-50 


0.1 


104 


-9.4199752 


244 


0.5 


-2.35 


1.0 


No 


50-50 


0.1 


133 


-13.899182 


245 


0.1 


-2.7 


1.0 


No 


50-50 


0.1 


1 


0.0000000 


246 


0.2 


-2.7 


1.0 


No 


50-50 


0.1 


25 


-2.0332112 


247 


0.3 


-2.7 


1.0 


No 


50-50 


0.1 


45 


-3.8604629 


248 


0.5 


-2.7 


1.0 


No 


50-50 


0.1 


61 


-5.6476472 


249 


0.1 


-2.35 


2.0 


No 


50-50 


0.1 


5 


-0.28666705 


250 


0.2 


-2.35 


2.0 


No 


50-50 


0.1 


44 


-3.6659021 


251 


0.3 


-2.35 


2.0 


No 


50-50 


0.1 


63 


-5.8169795 


252 


0.5 


-2.35 


2.0 


No 


50-50 


0.1 


62 


-5.7460627 


253 


0.1 


-2.7 


2.0 


No 


50-50 


0.1 


14 


-0.90617906 


254 


0.2 


-2.7 


2.0 


No 


50-50 


0.1 


22 


-1.6301615 


255 


0.3 


-2.7 


2.0 


No 


50-50 


0.1 


40 


-2.8798898 


256 


0.5 


-2.7 


2.0 


No 


50-50 


0.1 


38 


-2.7683023 


257 


0.1 


-2.35 


0.25 


No 


50-50 


0.1 


116 


-11.499357 


258 


0.2 


-2.35 


0.25 


No 


50-50 


0.1 


148 


-17.358444 


259 


0.3 


-2.35 


0.25 


No 


50-50 


0.1 


162 


-21.119200 


260 


0.5 


-2.35 


0.25 


No 


50-50 


0.1 


177 


-28.019414 


261 


0.1 


-2.7 


0.25 


No 


50-50 


0.1 


43 


-3.4548391 


262 


0.2 


-2.7 


0.25 


No 


50-50 


0.1 


87 


-7.2923632 


263 


0.3 


-2.7 


0.25 


No 


50-50 


0.1 


105 


-9.8528845 


264 


0.5 


-2.7 


0.25 


No 


50-50 


0.1 


135 


-14.206096 


265 


0.1 


-2.35 


1.0 


Yes 


50-50 


0.1 


41 


-3.1652049 


266 


0.2 


-2.35 


1.0 


Yes 


50-50 


0.1 


90 


-7.5318220 


267 


0.3 


-2.35 


1.0 


Yes 


50-50 


0.1 


114 


-11.360538 


268 


0.5 


-2.35 


1.0 


Yes 


50-50 


0.1 


146 


-16.893279 


269 


0.1 


-2.7 


1.0 


Yes 


50-50 


0.1 


3 


-0.11226917 


270 


0.2 


-2.7 


1.0 


Yes 


50-50 


0.1 


34 


-2.5620627 


271 


0.3 


-2.7 


1.0 


Yes 


50-50 


0.1 


53 


-4.4550489 


272 


0.5 


-2.7 


1.0 


Yes 


50-50 


0.1 


88 


-7.3062420 


273 


0.1 


-2.35 


2.0 


Yes 


50-50 


0.1 


6 


-0.28926393 


274 


0.2 


-2.35 


2.0 


Yes 


50-50 


0.1 


47 


-3.9974158 


275 


0.3 


-2.35 


2.0 


Yes 


50-50 


0.1 


67 


-6.2402162 


276 


0.5 


-2.35 


2.0 


Yes 


50-50 


0.1 


85 


-7.1175851 


277 


0.1 


-2.7 


2.0 


Yes 


50-50 


0.1 


8 


-0.46019542 


278 


0.2 


-2.7 


2.0 


Yes 


50-50 


0.1 


19 


-1.3618415 


279 


0.3 


-2.7 


2.0 


Yes 


50-50 


0.1 


33 


-2.5161509 


ZoU 


U.J 


-z. / 


z.u 


ICS 


JU-JU 


U. 1 


JO 




281 


0.1 


-2.35 


0.25 


Yes 


50-50 


0.1 


124 


-12.491142 


282 


0.2 


-2.35 


0.25 


Yes 


50-50 


0.1 


157 


-19.571515 


283 


0.3 


-2.35 


0.25 


Yes 


50-50 


0.1 


171 


-25.070431 


284 


0.5 


-2.35 


0.25 


Yes 


50-50 


0.1 


185 


-34.930124 


285 


0.1 


-2.7 


0.25 


Yes 


50-50 


0.1 


46 


-3.8622306 


286 


0.2 


-2.7 


0.25 


Yes 


50-50 


0.1 


97 


-8.5070409 


287 


0.3 


-2.7 


0.25 


Yes 


50-50 


0.1 


125 


-12.678477 


288 


0.5 


-2.7 


0.25 


Yes 


50-50 


0.1 


155 


-18.965882 



CE efficiency parameter 

b Exponent of the high-mass power law component of the IMF: |Kroupall2u0lT) (-2.35) or lKroupa & WeidneJ (20031) (-2.7). 

C Stellar wind strength parameter with which the "standard" ( Belczynski et al. 2010) stellar wind prescription is multiplied. 

^ Yes: all possible outcomes of a CE event with a HG donor are allowed, No: A CE with a HG donor star will always result to a merger. 

e Binary mass ratio distribution. "50-50" indicates half of the binaries originate from a "twin binary" distribution and half from flat mass ratio distribution. 

^ Parameter with which the "standard" iHobbs et al .1^20051) kick distribution is multiplied for BHs formed though a SN explosion with negligible ejected mass. 

° Ranking of the model in our statistical comparison. 

k Ratio of the likelihood of the observations given a model to our maximum likelihood model (logL(0|M, ) — log L(0\M re f)). 



